The pseudogap and photoemission spectra in the attractive Hubbard model 
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Angle-resolved photoemission spectra are calculated microscopically for the two-dimensional attrac- 
tive Hubbard model. A system of self-consistent T-matrix equations are solved numerically in the 
real-time domain. The single-particle spectral function has a two-peak structure resulting from the 
presense of bound states. The spectral function is suppressed at the chemical potential, leading to a 
pseudogap- like behavior. At high temperatures and densities the pseudogap diminishes and finally 
disappears; these findings are similar to experimental observations for the cuprates. 
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Real-space pairing [JD is the simplest physical idea that 
enables one to explain the pseudogap phenomenon ob- 
served for the normal state of high-T c superconductors 
(HTSC) PHl^l- At low temperatures and densities carri- 
ers are paired in weakly overlapping bound states sep- 
arated from the single-particle band by a binding en- 
ergy of the order of a few hundred degrees. The chem- 
ical potential is located between the two bands thereby 
reducing the intensity of the low-energy single-particle 
(photoemission, specific heat, tunneling), spin (suscepti- 
bility, nuclear relaxation rate), and particle- hole (optical 
conductivity) processes. At temperatures of the order 
of the binding energy, carriers become unbound and re- 
store the Fermi-liquid behavior. The layered structure 
of the cuprates supports the scenario, since reduced di- 
mensionality favors pairing. Also, the phenomenology of 
a charged Bose gas can be quite successfully used in the 
explanation of a number of normal and superconducting 
properties of HTSC 000]. 

It is therefore important to study model systems with 
pairing, such as the two-dimensional attractive Hubbard 
model. The simplicity of the model allows one to sepa- 
rate the net effect of the attractive interaction from the 
complications related to the origin of the pairing mecha- 
nism and to the complicated dependence of the effective 
potential on microscopic parameters. One successful ap- 
plication of the attractive Hubbard model to the physics 
of cuprates is due to Randeria and co-workers Jlml7|. 
Using the Quantum Monte Carlo method, they found 
a significant reduction of static spin susceptibility and 
nuclear relaxation rate at low temperatures and interme- 
diate couplings. Recently, Vilk et al. [M found a pseu- 
dogap in the spectral function of the attractive Hubbard 
model using Monte Carlo simulations and the maximum- 
extropy technique. Thus, it was demonstrated that real- 
space pairing can account for some unusual properties of 
HTSC. 

The self-consistent T-matrix approximation 19-21]] 



provides another method for studying dynamic proper- 
ties of the attractive Hubbard model. This approach is 
based on the low-density approximation to fermion sys- 
tems due to Galitskii [ p2| , which becomes exact in the 
limit of vanishing density. However, the resulting system 
of self-consistent integral equations is not easy to solve. 
Analytical treatment is very difficult, although it is some- 
times attempted [^3[|4). The full numerical solution of 
the equations is required. In the previous numerical stud- 
ies of the problem |2^,|2^] the equations were solved in 
imaginary times and results were then continued numer- 
ically to real times. In these papers relatively small cou- 
plings were studied and pseudogap features were found 
only at large momenta. 

In our previous paper p7| we formulated and 
solved numerically the T-matrix equations for the two- 
dimensional attractive Hubbard model for real times, 
thereby avoiding the necessity for analytical continua- 
tion. There, we focused on two-particle properties — pri- 
marily on the binding energy of pairs and its dependence 
on the particle density. In this paper we would like to 
discuss single-particle dynamics, in particular the single- 
particle spectral function which is directly related to 
angle-resolved photoemission spectra (ARPES). We find 
a clear pseudogap behavior of the ARPES at small mo- 
menta k, low densities n, and low temperatures T. With 
increasing n and T, the pseudogap disappears, in accor- 
dance with experimental observations for the cuprates 



The two-dimensional attractive Hubbard model is de- 
fined by the Hamiltonian 
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written in standard notation. e k — — 2i (cosfc^ + cosk y ) 
is the bare single-particle spectrum, \U\ is the coupling 
strength, and N is the total number of sites in the system. 
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The chemical potential fi determines the average particle 
density n. We regard Eq. (|l|) as a phenomenological 
model for the low-density system of holes in the normal 
state of HTSC. 

In the low-density limit n <C 1, one can make use of 
the small gas parameter and select only ladder diagrams 
in a diagrammatic representation of the T- matrix pO| , pl| , 
which leads to the expression 



T(q,w) = 
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where A(k, w) is the single-particle spectral function and 
/3 = (fc^T) -1 is the inverse absolute temperature. The 
real and imaginary parts of the self-energy S' and S" are 
expressed via T = T' + iT" as follows: 

E'(k lW ) = ^ / ^(q-k, Wl )x 

dw 2 / B (w 2 )T"(q,w 2 ) 



/F(wi)T'(q,w + wi) + 
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where /f,b(w) = [exp(/3w) ± are Fermi- and Bosc- 
functions respectively. Finally, the self-energy deter- 
mines the spectral function as 



A(k,to) 
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The integrals with singular kernels in Eqs.(^|) and ([|) are 
understood in the principal- value sense. The set of equa- 
tions (||)-@ is to be solved self-consistently for given 
|f7|, fi, and temperature T; then the particle density 
is given by n = 2N~ l Y, k J ^(k,w)/ F (w). Usually, 
a self-consistent solution is obtained iteratively, starting 
from a guessed form of A(k, to) and using the fast-Fourier- 
transform algorithm to calculate momentum-frequency 
sums p5| , p6| |. In our calculations we used a 64 x 64 lat- 
tice and a uniform mesh of 512 points in the frequency 
interval —20 1 < to < 30 1. The convergence of the itera- 
tive process is the major problem of the method, which 
puts limitations on the values of model parameters for 



which a self-consistent solution can be obtained. The 
convergence deteriorates for large \U\ and n and low T. 
The physically interesting values of \U\ start at ~ Qt, 
when the binding energy of the pairs is of the order of t. 
In this work, \U\ — 8t is used. For this coupling, itera- 
tions converge down to T = 0.3 1 for very low densities 
n < 0.03, and up to n ~ 0.20 for a high temperature 
T = 1.0*. 

Once a self-consistent solution is obtained, the inten- 
sity of the photoemission process is simply 



I(k,uj)=I {k)A(k,w)f F (u 
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where Io(k) involves the electron-photon matrix element, 
and is frequency- independent. Eq. (0) is approximate, 
for a discussion of its validity see, e.g., [^9|. In the fol- 
lowing we set io(k) = 1. 

In analysing the numerical results to be presented be- 
low, it is useful to keep in mind the exactly solvable 
atomic limit (t = 0) of the Hubbard model: 

i-A(k, to) = ^S(to + M +|C/|)+(l-|) S(to + M ), (8) 

from which the following properties are inferred, (i) The 
spectral function has the form of two peaks with weights 
§ and 1 — f (which are very different if n <C 1). (ii) 
The two peaks are separated by the binding energy of the 
pairs (which is \U\ in the atomic limit), (iii) At zero tem- 
perature, fj, = — 4^-, and in Eq. (0) the Fermi- function 
eliminates the second peak of A(k,u>). The resulting 
intensity J(k, to) is a single peak located ^ below the 
chemical potential. The system would therefore display 
a "pseudogap" (a true gap in this case) of size (iv) 
At higher temperatures, the Fermi-function is smoothed 
out giving rise to the second peak at a higher energy 
and weakening the first one, so the former might become 
stronger than the latter. 

In the general case of finite t/\U\, non-zero kinetic en- 
ergy leads to a number of new effects. It reduces the 
binding energy, i.e. the interpeak distance, and assigns 
finite widths to the peaks of A(k, to). Next, it restores the 
k-dependence of A(k, to) and I(k, to). Finally, due to the 
finite radii of pairs and their ovelapping, the binding en- 
ergy becomes density-dependent p^j . However, our nu- 
merical results show that properties (i)-(iv) listed above 
remain valid even at finite i/|E/|. Moreover, we believe 
they are generic to any fermion model with attraction in 
the low-density limit. 

In Fig. [I] we show the solution of Eqs. (0)-(R) for 
the lowest density n = 0.017 and T = 0.3 1. I(k,to) 
displays a complicated k- and w-dependence which can 
be understood as follows. A particle with momentum 
k can be found in two distinctly different states: either 
in a state of the single-particle band with an extended 
wave function or as a component of a bound state with a 
localized wave function. The two possibilities give rise to 
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FIG. 1. The intensity of ARPES I(k,w) = A{k,Lu)f F {u) 
for density n = 0.017, temperature T = 0.3 1, and several 
momenta k = m(j, j), m — 0,1,2, and 3 (from the top 
curve down). 

A(k, u>), which consists of two peaks separated by the pair 
binding energy AE = 2.1 t (for \U\ = 8t). Multiplication 
by the Fermi-function cuts off the high-energy (single- 
particle band) peak, which depends on the energy of the 
latter and the temperature. For k = (0,0) in Fig. |l| 
the high-energy peak is reduced in height significantly — 
to that of the low-energy peak, but not to zero. Note 
that after the cut-off, the peak is slightly shifted from its 
original position. The position of the high-energy peak 
disperses with k as does the bare specrum £k, and, as 
k increases, the peak gets cut off by the Fermi-function 
very rapidly (compare the cases for the different momenta 
k in Fig. |l| for u> > 0). Let us now turn to the low- 
energy peak. The probability of finding a particle with 
momentum k in a bound state is the square of the bound 
state's wave function. For zero total momentum P, one 
has 



V>p=(o,o)(k) = 



C 



E-2s(k)' 



(9) 



where C is the normalization constant and E the energy 
of the bound state measured from the bare atomic level. 
The relative height of low-energy peaks in Fig. [l] is in 
good agreement with iV^k)! 2 for E = 2e(0,0) — AE = 
— 10.lt. This corroborates the bound-state origin of the 
low energy peaks in A(k, to) and J(k, ui). Thus, on the 
basis of Fig. [|, we conclude that ARPES of the attractive 
Hubbard model exhibit a clear pseudogap behaviour at 
low temperatures and densities. The momentum and fre- 
quency dependences of the spectra have simple physical 
explanations. 

The temperature dependence of 7(k, w) for n = 0.017 
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FIG. 2. The temperature dependence of I(k,u>) for den- 
sity n — 0.017 and momentum k = (0,0). 

and k = (0,0) is shown in Fig. g. There are two main 
effects as temperature increases. First, weight is trans- 
ferred from the low-energy peak to the high-energy one. 
This is due to the progressive thermal excitation of par- 
ticles to the single-particle band and consequent weaker 
influence of the bound states on the single-particle spec- 
tral function. Secondly, the whole structure moves to 
higher energies relative to the chemical potential. These 
two effects lead to the rapid suppression of the pseudo- 
gap as temperature increases. Note that the distance 
between the two peaks is T-independent and remains ap- 
proximately the pair binding energy (slightly reduced by 
the cut off), in accordance with the atomic limit. 

Fig. H presents the density dependence of 7(k, ui) for 
T = 0.5 i and k = (0,0). Clearly, the pseudogap disap- 
pears as n increases. We have already argued elsewhere 
p7[ that this is a result of the rise of the two-particle 
level due to the packing effect when pairs begin to over- 
lap. (Intuitive arguments of this kind were given earlier 
in p!q|.) Since the binding energy decreases with n, the 
temperature becomes progressively more effective in un- 
bounding pairs, washing away the pseudogap. 

It is quite remarkable that such a simple system as 
the attractive Hubbard model and such complex systems 
as high-T c superconductors have very similar dynamical 
properties. They both display pseudogaps at low temper- 
atures and carrier densities, which disappear as T and n 
increase. This suggests the conclusion that the carriers in 
HTSC do experience some sort of short-range attraction. 
HTSC therefore exhibit properties which are generic to 
fermionic systems with attraction and which are captured 
in our model calculations. 

One could now proceed in elaborating the model while 
trying to keep the properties obtained intact. Further 
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insight into the problem can be gained by considering 
the opposite limit of nearly complete filling 2 — n <C 1. In 
this case, Eq. (]]]) may be viewed as a phenomenological 
model for electrons rather than holes. (To some extent, 
the nearly fully filled band imitates the nearly filled lower 
Hubbard band when the Coulomb repulsion is taken into 
account. Unfortunately, this analogy is not complete, 
due to different temperature behaviour of /i; see below.) 
The quantity J(k, w) of Eq. (Q) has now the meaning 
of the number of electrons emitted from the system per 
time unit, which brings the whole model closer to reality. 
There is no need to recalculate the spectra, since on a 
bipartite lattice the dilute and nearly filled limits are 
related by the particle-hole transformation, which loads 
to the relation 

A(k, w; n)f F (u) = A(k + Q, -w; 2 - n)[l - f F (-u)], 

(10) 

where Q = (it, ir) for the square lattice. Due to the in- 
version of the occupation numbers, it is now the pairing- 
induced peak that gets cut off by Fermi-function. The 
inversion of the frequency places the remaining peak be- 
low the chemical potential. The resulting ARPES are 
shown in Fig. [l| They are complementary to the spec- 
tra of Fig. |l|. Clearly, pseudogap is present, because the 
spectrum with the largest momentum k = (jr, it) is still 
peaked far below the chemical potential. 

The overall picture looks very much like spectra of a 
weakly-interacting system but shifted from fi by half of 
the pair binding energy. We emphasize that pairs them- 
selves are not seen explicitly in the spectra, since the pair- 
induced peak of A(k, uj) has been cut off by the Fermi- 
function. Nevertheless, the pairs are present implicitly, 



manifesting themselves in the shift of the chemical po- 
tential. This observation is important for understanding 
the ARPES of HTSC. 

We do not present temperature and density depen- 
dences of /(k, us) for the nearly filled case, for they are 
complementary to Figs. || and |j| respectively. The pseu- 
dogap now vanishes with decreasing electron density (in- 
creasing hole doping), in accordance with Fig. ^. The 
temperature dependence is, however, different from the 
dilute limit. At nearly full filling, the chemical poten- 
tial goes up, with temperature and the distance between 
(i and the single-particle band increases. Therefore, the 
pseudogap is expected to rise with T in this case. To ob- 
tain the correct temperature behaviour one would need 
to consider the density regime close to half-filling, which 
is outside the range of validity of the T-matrix approxi- 
mation. 

In conclusion, we have shown that photoemission spec- 
tra of the attractive Hubbard model in the low-density 
limit display a clear pseudogap behavior, qualitatively 
similar to that of high-T c superconductors. Our findings 
support the suggestion that the pseudogap feature ob- 
served in HTSC results from real-space pairing and the 
formation of bound pairs in the normal state. 
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subject. PEK acknowledges the support of EPSRC grant 
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